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We set up a simple transfer matrix formalism to study the existence of bound states at interfaces 
and in junctions between antiferromagnets and d-wave superconductors. The well-studied zero 
energy mode at the {110} interface between an insulator and a d x 2_ y 2 wave superconductor is spin 
split when the insulator is an antiferromagnet. This has as a consequence that any competing 
interface induced superconducting order parameter that breaks the time reversal symmetry needs 
to exceed a critical value before a charge current is induced along the interface. 

PACS numbers: 74.50+r, 74.72.-h, 74.25.Ha 



The discovery of the symmetry of the superconduct- 
ing order parameter has been one of the most successful 
studies of the High-T c materials. Angular resolved pho- 
toemission spectroscopy has revealed the nodes in the 
gap function and tunneling experiments have proven the 
sign, change between adjacent lobes of the d x 2_ y 2 wave 
gapErQ. It was first shown by Huo that this sign change 
can lead to zero energy Andreev bound states (ZEBS) 
at the surface of an insulator and a d-wave superconduc- 
tor. These Andreev bound states were later identified 
with the zero bias conductance peaks observed in tunnek 
ing experiments. The experiments by Covington et.ai.S 
indicated, however, that the surface states were sponta- 
neously split by a minigap. Several ideas were proposed 
for this effectEl; one of which included the induction of 
a time reversal symmetry breaking_js component of the 
order parameter near the interfaced. The resulting gap 
d + is lowers the condensation energy by lifting the di- 
rectional degeneracy of the ZEBSEj. Later Honerkamp 
et.alB used a tight-binding model with onsite repulsion 
and spin dependent nearest neighbor interaction to self- 
consistently study the competition between additional in- 
duced orders near the surface of an insulator and a d x 2_ y 2 
wave superconductor. 

Motivation for studying close domains of antifcrromag- 
netism and superconductivity arises from the existence 
of striped domains in the cuprate materials. This was 
further emphasized by recent elastic neutron scattering 
experiments showing that static antiferromagnetic order 
is induced in a superstructure around the vortices 4*1 
the mixed state of La2- a; Sr a ;Cu04li3 and La2CuC>4 +( 5liII. 
These experiments are consistent with a static environ- 
ment of alternating antiferromagnetic and d-wave super- 
conducting stripes around the vortex cores. Thus the 
electronic states in such an environment is an important 
question. 

Inspired by these experiments we set up a simple trans- 
fer matrix method to identify bound states on interfaces 
and junctions between antiferromagnets and d-wave su- 
perconductors. In particular we discuss a single interface 
separating antiferromagnetic and d-wave superconduct- 



ing half-planes (AF/dSC), and point out a few differences 
from the conventional non-magnetic insulator-d-wave su- 
perconductor interface (I/dSC). Note that the antiferro- 
magnetism forces us to study a lattice model which is 
contrary to the usual discussion of Andreev interference 
in terms of semi-classical continuum models. 
A simple lattice model that includes both d-wave super- 
conductivity and antiferromagnetism is given by the fol- 
lowing Hamiltonian 

H = -t ^ C Lcm<r + H.C. - fl^2 C l* c n* (1) 
(n,m)a no 

+ J2 A ^ c U c li + H - c - ( 2 ) 

(n,m) 

+ M " ( c 1t c «T - c Ii c «i) ( 3 ) 
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where (n, m) denotes nearest neighbors. M n and A„ jm 
are the spatially dependent magnetic and superconduct- 
ing order parameters. This Hamiltonian is quadratic and 
can be diagonalized by a Bogoliubov-de Gennes (BdG) 
transformation 

li = ^2 M") c Lx + o-w CT (n)c„_ CT . (4) 

n 

with <7 equal to +1 (-1) for spin up (down). We use the 
notational convention that the spin indices on u a and v a 
follow that on the Bogoliubov operators 7J. 
In the case of a <i a; 2_j / 2-wave superconductor there are 
two qualitatively different orientations of the interface; 
the {100} and {110} directions corresponding to a verti- 
cal and diagonal stripe respectively. Both cases are stud- 
ied below with the :r-axis (j/-axis) chosen perpendicular 
(parallel) to the interface which is placed at x = 0. The 
lattice constant is set to unity. Assuming translational 
invariance along the y-direction the AF/dSC interface 
reduces to a one dimensional problem. For the {100} 
interface the resulting Bogoliubov-de Gennes equations 
have the form 
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e a u qa (x) = -t (u qa (x + 1) + u qa (x - 1) + 2cos(q)u qa (x)) - fiu qa (x) + aM x u q+Qa (x) 
+ i^i+i.x) v Q * (x + 1) + (Af _ liX ) v qa {x - 1) + 2 cos(g) (-A£) w 9tr (x) 

finger (x) = i (w gcr (x + 1) + u gCT (x - 1) + 2 cos(<7)v 9(T (x)) + jtii; g(r (x) + aM x v q+Qa (x) 
+ (A*i 1;X ) <v (» + !) + (^-i,*) *V (a; - 1) + 2 008(9) {-K d ) *V (x) 



(5) 
(6) 



after fourier transforming along the y direction. The cor- 
responding equations for the fourier components u q +Q a 
and v q +Q a are obtained by simply performing the substi- 
tution q — > q + Q. These BdG equations are diagonal in 
the spin index with the only difference between spin up 
and down being the sign of the magnetic term. 
A simple way to study bound states aLthe interface is in 
terms of the transfer matrix methodE3. Thus we intro- 
duce a (q, e)-dependent matrix T (x + 1, x) defined by 

*(x+l) =T(x+l,x)*(x). (7) 

which transfers the spinor W from site x to site x + 1. 
For a model with nearest neighbor coupling ^> takes the 
explicit form '5 (x) = (ip (x) ,tJj (x — 1)) where 

1p (x) = (u qa (x) , V qa (x) , U q+ Q„ (x) , V q+ Qo (x)) . (8) 

The associated 8x8 transfer matrix has the general form 

T(x+l,x) = (f^ (9) 

where A (B) denotes the 4x4 coefficient-matrix con- 
necting -0 (x + 1) and tp (x) (jp (x — 1)) determined from 
the BdG equations (|[]|) . In the simplest case of a sharp 
interface we have the following spatial dependence of M x 
and A x 

M x = M(-l) x for x<0 (10) 
A x = A d for x > (11) 

Thus there are effectively three different transfer matri- 
ces; one in the bulk magnetic region Tm, one in the 
bulk superconducting region T$c an d one associated with 
transfer through the interface Tj. By diagonalizing Tm 
and Tsc there exists decaying, growing or propagating 
eigenstates depending on whether the eigenvalues are 
less, larger or equal to one, respectively. Here, decaying 
and growing are referred to propagation along the x-axis 
for increasing x. If PETm denotes the matrix obtained 
after propagating the eigenvectors of the bulk magnetic 
transfer matrix through the interface we introduce a ma- 
trix a given by 

PET M = ETsc-a (12) 

where ETsc is the matrix containing the eigenvectors of 
the bulk superconducting region as coloum vectors. The 



dot indicates matrix multiplication. Now, let S 1 ™ and 
Sg C denote the subspace of growing eigenstates of PETm 
and ETsc respectively, and consider the following linear 
combination of the growing states of PETm 

]T fc\PET M i> = J2J2 p i a j i\ET sc j> (13) 

= E E \ET SCJ > 

From equation ([l3]) it is evident that to have a bound 
state at the interface the vector j3 must belong to the null 
space of the reduced matrix a r , which is the Sg C x S 7 g n 
upper left part of the original matrix a since the matrices 
PETm and ETsc are organized to have the eigenstates 
with the largest eigenvalues as coloum vectors to the left. 
In the case that the two subspaces Sg C and 5™ have the 
same dimension a bound state at the interface is charac- 
terized by the vanishing of the determinant of a r 

Bound states: det (a r ) = 0. (14) 

Plots of the wavefunctions with values of (q, e) that sat- 
isfy Eqn. (|l4|) verifies that these states indeed are bound 
to the interface (not shown). The following explicit val- 
ues of the input parameters are chosen: t = 1, A<j = 0.14, 
M = 2.0 and /i = —0.99 (for simplicity we ignore next- 
nearest neighbor coupling). Figure la shows the determi- 
nant plottet as a function of energy for the {100} inter- 
face. There are bound states close to the superconduct- 
ing gap edge that disperses downward in a cosine form 
(Figure lb). These are the well-known de Gennes/Saint- 
James states existing, on the surface of an insulator and 
a superconductoilijo. 
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FIG. 1. {100} interface between an antiferromagnet and a 
d-wave superconductor: a) Determinant of a r as a function 
of energy e for q — 0.1. There is a de Gennes/ Saint- James 
bound state close to the superconducting gap edge. As seen 
in b) , their dispersion has the expected downward cosine form 
until it merges with the continuum. 
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The induction of additional gap symmetries, extended 
s-wave or p-wave, near the {100} interface of a d-wave 
superconductor and an antiferromagnet has been studied 
self-consistently by Kubokilj. These local gap perturba- 
tions will slightly modify the graphs in Figure 1. There 



is no spin splitting of the dGSJ mode in this geometry. 
We turn now to the more interesting configuration of 
a {110} interface. Allowing for a possible interface in- 
duced sub-gap order with extended s-wave symmetry the 
Bogoliubov-de Gennes equations have the form 



i 

e a u qa (x) = -2f cos(p) (u qa (x + 1) + u qa (x — 1)) - [iu qa (x) + oM x u qa (x) (15) 

- 2i sm{q) (A d x+hx v qa (x + 1)- A d x v qa (x - 1)) + 2^ cos(g) (A x+1>x v qa (x + l) + A s x v qa (x - 1)) 
^<yV qcr (x) = 2t cos(p) (v qa (x + 1) + v qa (x — 1)) + liv qa (x) + oM x v qa (x) (16) 

- 2i sin(g) (A* x d +lx u qry (x + 1) - A x d u qa (x - 1)) - 2^ cos(g) (A* x % lx u qa (x + 1) + A x s u q<J (x - 1)) 



These equations are diagonal in the fourier component q 
obtained after fourier transforming parallel to the {110} 
interface since there is no staggering of the moments 
along a diagonal line in a square antiferromagnetic lat- 
tice. In Figure 2 we plot again the determinant of the 
reduced matrix a r as a function of energy e when A s = 0. 
As seen the spin degeneracy of the ZEBS (dashed curve) 
is lifted at a {110} AF/dSC interface. As opposed to 
the usual dGSJ states in Figure 1, this splitting is also 
caused by the fact that a {110} interface belongs to only 
one sublattice whereas the {100} interface studied above 
contains the same amount of spin up and down sites. 



acts as a local magnetic field. A similar effect caused by 
a correlation induced magnetization near the interface in 
the case of a I/dSC surface was discussed by Honerkamp 
et.aB. This "Zeeman" effect is also directly related to 
the split zero energy Andreev mode observed in the cen- 
ter of vortex cores of underdoped cuprateSj-whnre local 
antiferromagnetism has been shown to existt3~E3. 
To the best of our knowledge there has been no self- 
consistent calculation investigating any {110} AF/dSC 
interface induced subdominant order parameters. 
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FIG. 2. Determinant of a r versus the energy e for the {110} 
AF/dSC interface. Again this is plottet inside the supercon- 
ducting gap and with q — 0.1. The dashed curve is the usual 
case of an I/dSC interface which clearly contains a ZEBS 
(the insulator state is obtained by performing the substitu- 
tion M n — * —M n for the hole part of the BdG equations 
only). The solid curves show the spin splitting of the ZEBS 
for this particular value of q. 

The splitting of the ZEBS by A s -mixing in the usual 
situation of a I/dSG interface has been extensively stud- 
ied in the literatureQ~El It is also well-known that a mag- 
netic field further splits the ZEBSI3. The above spin split- 
ting at AF/dSC interfaces is similar to this magnetic field 
effect in the sense that the magnetic interface effectively 
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FIG. 3. a) Same as in Figure 2, but with an induced ex- 
tended s-wave gap function near the interface, i.e. d — » d + is. 
For clarity we do not show the original ZEBS (dashed curve 
from Fig. 2). b) Schematic representation of the splitting of 
the original zero energy Andreev bound state (dashed curve): 
1) The antiferromagnetic interface breaks the spin degeneracy, 
as shown in Fig. 2. 2) Induction of a possible sub-dominant 
s-wave gap parameter A s further splits the spin up/down 
states by breaking the directional degeneracy. 3) Only when 
A s exceeds a critical value is an interface current induced. In 
this last figure, which corresponds to the situation from a), 
A s is equal to A d on the interface and decreases linearly to 
zero within 20 sites of the interface. 

However, we know from the study of I/dSC surfacesi'0 
that the strong pair breaking effects of a {110} geometry, 
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as opposed to a {100} surface, tends to stabilize the sub- 
dominant s-wave component. Thus, even though there 
is no Fermi surface instability begging for removal of the 
ZEBS from the Fermi level in the case of a AF/dSC {110} 
interface, one should still consider the effect of an addi- 
tional local superconducting order parameter is compet- 
ing with the splitting caused by the magnetism. The 
consequences of this competition for the ZEBS are dis- 
cussed in Figure 3. 

The induction of a surface current is a well-known con- 
sequence of the tffue reversal symmetry broken state of 
I/dSC interfacesQft However, for the AF/dSC interface 
with a locally induced d ± is order parameter there is 
a criticaJ-jvalue of A* before a current runs along the 
interfaces. In Figure 3a we show the situation when the 
induced A s has exceeded this critical value. Figure 3b is 
a schematic representation of the splitting of the original 
ZEBS with the first sketch corresponding to the parame- 
ters from Fig. 2 and the last sketch to those from Fig 3a. 
We stress that only a self-consistent model calculation 
can determine the magnitude of the directional splitting 
caused by is compared to the spin splitting caused by 
the antiferromagnetism, and hence the relevancy of the 
interface current. 

In conclusion we have set up a simple method so deter- 
mine the existence of bound states at the interfaces of 
d-wave superconductors and antiferromagnets. In par- 
ticular we studied the energetics of the notorious zero 
energy mode Ippund to {110} I/dSC interfaces first dis- 
covered by HvJj. This state is always spin split when the 
insulator is an antiferromagnet and is analogous to the 
split states found around the magnetic vortex cores of 
YBCO and BSCCO crystals. In the case of an array of 
junctions corresponding to a periodic domain of vertical 
or diagonal stripes these states will hybridize and even- 
tually form a band. A current along the interface exists 
only when the effect of a competing, interface induced is 
component exceeds the spin splitting. 
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